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ABSTRACT 

It has been suggested that the anomalous Na, Mg and Al observed in Globular 
Cluster Red Giant stars could be the result of a thermally unstable hydrogen 
shell. Currently accepted reaction rates indicate that temperatures of approx- 
imately 70-75 million K are required to produce the observed enhancements 
in Na and Al along with depletions in Mg. 

The work presented here attempts to model the H shell instability by a 
simple mechanism of altering the energy production in the region of the H 
shell. We show that even extreme cases only give rise to small intermittent 
temperature increases that have minimal affect on the surface abundances. 
Full evolutionary modelling incorporating this technique simply accelerates 
the evolution of the RGB phase producing the same surface abundances as 
other models but at an earlier time. We conclude that, unless hydrogen shell 
instabilities manifest themselves quite differently, they are unlikely to lead 
to the required temperatures and alternative explanations of the abundance 
anomalies are more promising. 
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1 INTRODUCTION 

As long ago as 1947, Popper ( |1947| ) had noted that the globular cluster red giant L199 in 
M13 was CN-strong. Further work by Harding ( |1962| ), Osborn ( |1971| ) and Hesser, Hartwick 
and McClure ( |197(j| ) confirmed that this anomaly is common. So for over fifty years it has 
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been known that the surface abundances of elements such as C and N can vary from one 
red giant star to another within an individual globular cluster. It has become clear in the 
intervening years, however, that these anomalies are not restricted to C and N: heavier 
elements such as O, Na, Mg and Al also show star to star variations. Most other elements 
studied show little variation, and the degree of variation of C, N, O etc. differs from cluster 



to cluster and is much less apparent in field stars (|Langer et. al.l99^) . 

Observations by many groups, including Kraft et al. (|1997|) , Smith and Kraft (|1996D , 
Langer (|1992| ), Shetrone (|1996a|) ( |1996b|) (|1997| ), have identified this abundance anomaly 
problem and the most commonly suggested explanation for such high Na and Al are the deep 
mixing ( ^weigart fc Mengell979|) and primordial enhancement (|Cottrell fc Da Costal98TD 
hypotheses. The deep mixing hypothesis suggests that the rotation of the star can give rise 
to meridional circulation currents that in turn create a mixing zone in the radiative region 
separating the top of the hydrogen burning shell (HBS) and the base of the convective 
envelope (BCE). This would allow elements processed in the HBS to be mixed to the surface 
as the star proceeds up the RGB. The primordial enhancement hypothesis suggests that the 
stars are born with these anomalies already present. The stars are formed from material 
ejected by intermediate mass 3 — IOMq) AGB stars from an earlier epoch. Without 
going into the full mechanisms of each process, either could in principle explain the observed 
anomalies, but modelling and observations have indicated that probably both processes occur 
to some extent. 

Langer, Hoffman and Zaidins ( [1997|) suggested that the Na and Al abundance anomaly 
problem may be solved if the H shell reached higher temperatures for at least part of its 
lifetime whilst on the RGB. The idea stems from work by Von Rudloff and Vandenburg 
( |1988|) and even earlier work by Bolton and Eggleton ( |1973| ). Von Rudloff and Vandenburg 
( |1988|) found that the H shell of a non-rotating small mass red giant model was stable but 
only marginally. A rotating model could be unstable. The link with rotation is particularly 
encouraging because this mechanism could explain both the deep mixing and the thermally 
unstable shell. The work performed by Langer et al ( |1997| ) shows that a thermally unstable 
H shell would allow for the greater production of ^^Al at the expense of ^^Mg and that there 
would also be extra ^'^Na available for mixing to the surface. 

The conclusion was that if the H shell reached temperatures of approximately 70-75 
million K then the material mixed to the surface would produce a composition in close 
agreement with observation. In fact temperatures of this order and no lower would be re- 

© 1994 RAS, MNRAS 000, §-0 



Stellar abundance anomalies and H shell instabilities 3 

quired to match the Na observations because at the lower temperatures too much ^^Na 
would be produced. Also it is only at these temperatures that significant ^'^Al is produced 
and ^^Mg is only destroyed at temperatures above 70 million K. Therefore, to have any 
surface depletions in ^^Mg the temperature of the H shell has to reach at least 70 million K 
unless primordial changes are invoked or the accepted reaction rates are altered. The pre- 
ferred temperature range was 72-73 million K (see section 2 of Langer et. al ( p.997| )). This is 
the temperature at which they believed all the observed surface abundance anomalies could 
be matched. 

As the structure of a model is heavily dependent upon the temperature profile any change 
to the temperature must be included in the evolution modelling. The following sections 
discuss models where temperatures in the H shell are artificially increased by decreasing 
the energy generation rate. The next section looks at oscillating temperature rises in the 
HBS; i.e. a pulsating model where each pulse immediately follows the previous pulse. In the 
following section pulses in the temperature of the HBS are created with a time gap between 
each pulse. 

2 COMPUTER CODES 

The calculations reported here make use of two separate computer codes: the Monash/Mt Stromlo 
stellar evolution code, with which we model the evolution of stars from contraction to the 
main sequence up to the RGB tip, and a post-processing nucleosynthesis code ( |Cannonl990D , 
with which we model the abundances within the star for the same epoch. The nucleosynthe- 
sis code uses the independent variables from the evolution code for its structural basis, with 
typically every sixth evolution model being used. The full details of how these two codes 
work and interface can be found in Messenger ( pOO0|) . 

The postulated deep mixing has been included in the nucleosynthesis code, and is de- 
scribed by three parameters - the depth, the speed and the time of onset of the deep mixing. 
We measure the depth in the same way as Wasserburg, Boothroyd and Sackmann (|1997D , 
by mixing down to a temperature AlogT above the base of the HBS, which is defined as 
the first point (moving outward from the center) to have a non-zero hydrogen abundance. 
This mixing commences once the HBS burns through the discontinuity left from the first 
dredge-up episode earlier in the evolution. Prior to this time it is assumed that this molec- 
ular weight discontinuity will inhibit mixing: once it is removed, the mixing is assumed to 

© 1994 RAS, MNRAS 000, 



4 B. B. Messenger and J. C. Lattanzio 

begin. Finally, we found (in agreement with previous researchers) that provided the speed 
of mixing was above a small critical value, it made almost no difference to the abundance 
patterns. We used a mixing speed of lO~^M0/year in all calculations reported here. 

In the work of Messenger ( |2000| ) it was found that varying the mixing depth as the star 
ascends the RGB gives better agreement with observations of Na. Therefore this type of 
modelling is adopted in this paper, and the mixing depth was changed linearly from an 
initial setting A logT = 0.18 to A logT = 0.03 as the star ascended the RGB. 

3 INCREASED H SHELL TEMPERATURES 
3.1 Algorithm 

To increase the temperature of the H shell it is not simply a matter of arbitrarily adding 
temperature to the mesh points within the shell. The temperature is one of the dependent 
structure variables and affects other variables such as the density and pressure through 
the equation of state. To self-consistently increase the temperature we chose to change the 
energy generation of the shell by increasing or decreasing the energy production rate for 
each point if it is within the boundaries of the H shell. 

An increase of the energy generation rate causes the star to expand and the temperature 
to decrease. A reduction in the energy generation rate causes the opposite; the star contracts 
and the temperature increases. Therefore it was necessary to decrease the energy rate from 
H burning in order to artificially create a temperature increase. 

A factor for the energy generation rate reduction, designated as was introduced. This 
reduction is only initiated after the H shell passes through the /U- discontinuity remaining 
from the first dredge-up episode. The desired nuclear energy generation rate at a point in 
the H shell (when this feature is switched on) is 

^used = erate(l — d) (l) 

However, it is not simply a matter of reducing the energy generation rate over the full 
shell. This would create a large jump in energy generaqtion from the point just below the 
shell and also from the point at the top of the shell.|^ It is necessary to linearly interpolate 

* We define the top of the shell to be the point that first has a hydrogen abundance 0.98 of the surface abundance moving 
from the centre of the model outwards. 
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Figure 1. Basic energy generation rate reduction algorithm. The mesh-based energy generation reduction factor hm is linearly 
increased from zero to g over n mesh points, held at g for 2n mesh points and then linearly decreased to zero over a further n 
mesh points. 



from to the value of g over a number of mesh points at the bottom of the shell and from 
(7 to at the top of the shell. 

Let hm be the mesh-based linear interpolation factor. This varies from zero at the bottom 
of the HBS to a maximum of g over n mesh-points, remains at this maximum over 2n mesh 
points, and then decreases linearly over the next n mesh points to a value of zero again. 
This is shown in Figure 0. 

Once hm has been determined the energy generation rate remaining at any point within 
the H shell (as long as the H shell instability feature is switched on) is 



This alone caused a temporal discontinuity in the models and subsequently created con- 
vergence problems. Remember that this arbitrary temperature rise is non-physical. To avoid 
these convergence problems it was necessary to smoothly introduce the reduced energy gen- 
eration rate in time as well as in mass (i.e. to gradually turn on the energy reduction). 
We introduce the reduction rate gradually over time with a parameter ht which is linearly 
increased from to 1 over 20 timesteps. Thus the energy generation rate within the H shell 
is 



Note that parametrisations in terms of mesh point (rather than mass) and model number 
(rather than time) are purely numerical conveniences. In reality, of course, these variations 




(2) 




(3) 
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Figure 2. Model for a constant energy generation rate reduction. As can be seen the greater the reduction of the nuclear 
burning energy rate, the more rapid the evolution of the model. For g = 0.5 the time spent on the RGB (since first dredge-up) 
is approximately 7.3% less than a star where 9 = 0. Similarly for g = 0.8 the time spent is approximately 12% less and for 
g = 1.0 the time spent is approximately 16% less. 

should be dependent on mass and time, but without any identified mechanism for the in- 
stabihty, and in view of the exploratory nature of the calculations reported here, we believe 
this simplification is justified and in no way affects our conclusions, which depend only on 
the shell temperature. 



3.2 Results 

All calculations reported here are carried out for a star of 0.8Mq with a metallicity of 
Z = 0.0005. Parameters relevant to the deep mixing were given above. 

The temperature in the H shell is determined by the requirement that the energy pro- 
duction rate be high enough to support the outer layers of the star. With a decreased rate of 
energy generation at a given temperature, the star's response is to increase the temperature 
in the H shell. As the star evolves up the RGB, the H shell progresses outwards and thins. 
The higher temperatures mimic this process and the final result is an acceleration of the 
star's evolution. This is shown in figure ^, where plots of luminosity against time for differing 
values of g are shown. The greater the energy rate reduction, the quicker the star evolves to 
the RGB tip. 

This result will not assist in the abundance anomaly problem because any enhance- 
ments/depletions will still occur to about the same degree as before but just at an earlier 
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epoch. The higher temperatures within the H shell mimic normal evolution but reach it at 
an earlier stage. All this manages to do is force the model to the He flash more quickly. No 
extra Na, Mg or Al would be produced in this situation, in fact, with a reduced timescale 
for the burning, the effect may be the opposite of what is required, but in any event it will 
be quantitatively negligible. 



4 INTERMITTENT TEMPERATURE RISES 

Langer et. al. ( p.997| ) suggested that the H shell is a thermally unstable region and that 
under certain conditions the shell may give rise to higher temperatures than standard physics 
suggests. The modelling performed in the previous section is not that of a thermally unstable 
region, but just a consistently higher temperature. To model a thermal instability the energy 
generation rate modification must not be steady but should be periodic or at least variable. A 
small modification to the algorithm of section ^ was therefore made to simulate a thermally 
unstable H shell. 

The simplest way to do this is to continue to use the time based rate reduction algorithm 
but just make it periodic. Here the energy generation rate reduction was introduced over five 
timesteps, held at the peak for 20 timesteps and then decreased over another five timesteps. 
This was then repeated until the RGB tip was reached. For an initial test and to determine 
how pulsation may affect convergence there was no gap between each pulse. Apart from the 
increase and decrease of each pulse the affect would be identical to the continuous algorithm 
described in the previous section. We would therefore expect little change from the model 
for the previous algorithm and that is exactly the case. Figure |^ which has g = 0.8 shows 
that the result of the above algorithm is to accelerate the evolution just as it did for the 
case of a continuous H shell temperature rise. The problem with the above approach is that 
the decrease in the energy generation rate causes the shell to contract and temperatures to 
rise. This shifts the star to a different evolutionary track. Once the energy generation rate 
returns to normal values the shell begins to expand again but there is not enough time for 
the star to settle back to its original (or close to original) evolutionary track before the next 
pulse begins. The end result is that the star shifts from its original evolutionary track and 
diverges from it henceforth. 

If we do not want to significantly alter the evolution of the star, then it is necessary to 
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Figure 3. A periodic energy generation rate reduction model and a standard model with normal energy production. 

introduce a time delay between pulses which is long enough for the star to resettle to its 
original track. 

We use a simple method, but it is efficient and can test the concept of a H shell instability 
without the introduction of large amounts of complex code which is hard to justify since 
there is no conffimed mechanism for the instability. The purpose of our algorithm is to test 
the principle of a periodic H shell instability, not the details. We use a technique similar in 
principle to that shown in figure 1 but we introduce periodicity and time delays between 
pulses. We ffistly select the desired number of pulses between FDU and the tip of the giant 
branch. We thus determine the inter-pulse duration and distribute the pulses evenly over the 
giant branch life-time. After some experimentation we decided to use 20 pulses because this 
allows the star to settle back to its standard evolutionary track and it should illustrate any 
differences likely to occur due to the pulses. We then take as a characteristic time-scale the 
time-step used to construct the current model, which we denote as dt. Next, the duration of 
a pulse is taken as 30(it, comprising bdt over which we linearly increase ht from zero to unity, 
then 20(it using ht = 1, and a further 5dt when hf is linearly decreased back to zero. Of 
course, the time-steps used by the code vary, and are not all equal to dt during the pulse, but 
dt gives a convenient timescale for the test. Also, because the time-steps vary with evolution 
along the giant branch, the dt used for each of the 20 pulses will be slightly different, but 
this is not important for our exploratory calculations. 
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The result of the introduction of this algorithm is shown in figure ^. The points that 
should be noted are the decrease in the luminosity shown during a pulse accompanied by 
a temperature increase and the fact that the curve and hence the model settles back to its 
original track. It is only at the latter stages of the RGB evolution that the curves begin to 
diverge. 

The conclusion drawn from this model is that the H shell temperature can be increased 
intermittently without significantly affecting the evolution of the star. But how will it affect 
the nucleosynthesis of the model? Is the temperature increase high enough and long enough 
to produce significant amounts of ^^Al and reduce ^^Mg as originally hoped? It is necessary 
to look at the temperatures within the H shell to see whether temperatures of approximately 
70 million K are reached, the temperature at which ^''Al enhancements are reasonable and 
^^Mg depletions begin. Figure |^ shows the temperature at the base of the H shell as the 
star ascends the RGB for an energy generation rate reduction of 50% {g = 0.5). The peak 
temperature reached by any of the pulses is only increased very slightly, well short of the 
required temperatures for ^"^Mg destruction. Even a 100% reduction of the nuclear burning 
energy generation rate {g = 1.0) during a pulse, as shown in figure |], cannot produce high 
enough temperatures within the shell for the processes discussed here. For a test case, some 
of the gravitational and neutrino energy generation rates were also reduced {g = 1.5). This 
created a minor increase in the temperature (see figure but the peak temperature still 
fell short of the required temperatures. 

Also if there are to be enhancements from this temperature increase, they must reach such 
temperatures early enough for significant composition differences to be seen in the envelope. 
The peak temperatures from these models are only seen during the very latest stages of 
RGB evolution, probably not early enough to match the data from M13 for example. 

It seems that the required temperature increase cannot be produced via this method. The 
suggested instability has some desirable features but to self consistently produce the required 
temperatures proves difficult: our simulated instability does not succeed, and the question 
of whether a genuine instability would succeed remains open until such a mechanism can 
be found, and consistently modelled. We are therefore reluctant to dismiss this possibility, 
based on the technique attempted in this paper. Also a genuine instability may be different 
to what we have modelled here and may create higher temperatures in the HBS. 

Despite this, temperature rises are seen whilst the pulse is active, and these could alter 
the nucleosynthesis. It would not be expected that any ^''Mg would be destroyed at these 
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Figure 4. Periodic energy generation rate reduction model and a standard model with no energy generation rate modifications. 
Here the evolution is not accelerated as much because there is enough time between pulses to allow the model to resettle on to 
its original evolutionary track. The model does deviate as it approaches the RGB tip but this could be avoided by selecting a 
smaller number of pulses. 
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Table 1. Reduced energy generation rate models. 

temperatures, but ^''Al production starts at lower temperatures and peak ^^Na production 
at even lower temperatures. Even though the temperature pulses are not high enough to 
test the hypothesis espoused by Langer et al (|1997|) , there may still be some effect on the 
surface abundances of Na, O and Al as well as of the other lighter elements including such 
isotopes of Mg as ^^Mg and ^^Mg. 

Table |l| lists all the models run where the energy generation rate was reduced. Case 6 was 
run through the nucleosynthesis code with deep mixing included. Figures pi p| and llO show 



the surface abundances for Na, Mg and Al respectively with the curves for the standard 
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Figure 5. The temperature at the base of the hydrogen shell for a periodic energy generation rate reduction model with 
9=0.5. 
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Figure 6. The temperature at the base of the hydrogen shelU for a periodic energy generation rate reduction model with 
9=1.0. 

model (Case 0) superimposed for comparison. It is quite clear that the temperature pulses 
have had a negligible effect on the surface abundances. 

The temperature rises are either too small or too short lived to have any effect. Longer 
lived or more pulses only accelerate the evolution of the star and does not assist in main- 
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Figure 7. The temperature at the base of the hydrogen shell for a periodic energy generation rate reduction model with 
3=1.5. (Thus removing all the burning energy, some, possibly all, the gravitational energy and the total energy may even be 
negative.) 
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Figure 8. [Na/Fe] versus [O/Fe] for Case 6 in Table 1 (solid curve) and a standard model with no energy generation rate 
modifications (dotted curve, Case 0). The data points are for M13. 
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Figure 9. Same as figure ^ but for [Mg/Fe] versus [O/Fe]. 




Figure 10. Same as figure ^ but for [Al/Fe] versus [O/Fe]. 
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Figure 11. [Na/Fc] versus [O/Fe] for Case 2 (solid curve) and a standard model with no energy rate modifications (Case 0, 
dotted curve). 

taining higher temperatures and thus providing more ^^Al for example. To test this, one of 
the models with constant energy generation rate reduction, Case 2 from table |l] was also run 
through the nucleosynthesis code, with its results and that of a standard model shown in 
figures |TT], [T^ and |TB|. Once again there is little effect on the surface abundances, as expected, 
although the Na values fit a little better. The process accelerates the evolution but does not 
significantly alter the surface abundances. 



5 CONCLUSION 

It appears that a thermally unstable H shell, modelled with the methods of this paper, 
cannot produce the high temperatures of about 70 million K required to produce the observed 
aluminium abundances. We succeeded only in accelerating the evolution with negligible effect 
on the interior temperatures and the surface abundances. Although a genuine instability, 
may manifest itself in a way which is quite different to these exploratory calculations, we 
have tried to artificially reproduce the most favourable conditions without success. 

It appears that primordial abundances variations among the stars combined with deep 
mixing currently offer greater potential in making progress on this problem. 
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Figure 13. Same as figure but for [Al/Fe] versus [O/Fe 
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